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Abstract 



The structural stability and electronic properties of the ternary inter- 

o 

^D \ metallic compounds NiSnM (M = Ti, Zr, Hf) and the closely related Heusler 



compounds Ni2SnM are discussed using the results of ah initio pseudopoten- 
tial total energy and band-structure calculations performed with a plane wave 
basis set using the conjugate gradients algorithm. The results characterize the 
lowest energy phase of NiSnM compounds, with a SnM rocksalt structure sub- 
lattice, as narrow gap semiconductors with indirect gaps near 0.5 eV, while 

><" 
Vh ' the Ni2SnM compounds are described as normal metals. Two other atomic 

C^ \ 

arrangements for NiSnM in the MgAgAs structure type result in energetically 

unfavorable compounds which are metallic. The gap formation in the lowest 
energy structure of NiSnZr and relative stability of the three atomic arrange- 
ments are investigated within a tight-binding framework and by considering 
the decompositions of each ternary compound into a binary substructure plus 
a third element sublattice. The stabilization of the lowest energy phase of 
NiSnZr is found to be mainly due to the relative stability of the SnZr rock- 
salt substructure, while the opening of the gap induced by the addition of 
the symmetry-breaking Ni sublattice makes a relatively minor contribution. 
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The results from the theoretical calculations for the NiSnM compounds are 

compared with the existing experimental data. From analysis of structural 

and chemical trends in the NiSnM compounds, CoVSn is predicted to be a 

new semiconducting intermetallic compound in the MgAgAs structure type. 

Preliminary first principles calculations suggest an indirect gap of 0.8 eV. 
61.50.Lt, 71.25.Pi, 71.25.Tn, 61.44.+p 
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I. INTRODUCTION 

Recent discoveries of unconventional electronic properties in multinary materials with 
complex crystal structures, such as quasicrystals and high T^ superconductors, have stimu- 
lated a growing demand for realistic theoretical studies of these systems. A full understand- 
ing of the interrelationships between composition, structure, and electronic properties of 
these complex systems requires a detailed analysis of individual materials at the microscopic 
level. While the quantum mechanical modeling of solids through first principles calculations 
has served as a very useful tool to obtain the microscopic information needed for struc- 
turally and chemically simple solids, there are practical problems associated with complex 
systems, especially with nonperiodic structures like quasicrystals. The challenge, therefore, 
is to identify simpler systems which can be fully studied with this theoretical tool, and yet 
shed light on issues that are relevant to the physics of more complex systems. 

One issue that arises in the study of complex intermetallic systems is the formation of 
a (pseudo)gap in the electronic density of states (DOS) at the Fermi level and the relation 
of this feature to the stability of the particular structure. For quasicrystals, the presence 
of an (almost) universal suppression in the electronic DOS at the Fermi level has long been 
established both experiment allytJ and from first principles band-structure calculations for 
rational approximants.a It is widely believed that this particular feature of quasicrystals is 
related to their stabilityB based on an analogy to the Hume Rothery rulea for the stability 
of metallic alloys. The role of a pseudogap in structural stability has also been suggested 
by the experimental observation of very high resistivities in structurally well ordered stable 
quasicrystals.™ Motivated by these observations, there has been recent interest in under- 
standing the structural energetics, stability and transport properties of intermetallic com- 
pounds with small unit cells exhibiting a (pseudo)gap in their electronic spectrum. However, 
relatively few examples of such systems are known. CsAu is an extensively studied exampldll 
with a band gap of 2.6 eV. The electronic structure and transport properties of AI2RU have 
recently been of considerable interest due to its structural simplicity and compositional 



similarity to quasicrystals.&u Another interesting class of intermetallic compounds with a 
pseudogap at the Fermi level is the skutterudite pnictides MA3 (M = Co, Rh or Ir, and A 
= P, As or Sb) . These compounds exhibit unusual transport properties making them likely 
candidates for advanced thermoelectric materials.^ The identification of additional examples 
of intermetallic compounds with (pseudo)gaps is therefore a problem of both fundamental 
and technological interest. 

A few years ago, the ternary intermetallic compounds NiSnM (M = Ti, Zr, Hf) were re- 
ported to exhibit unusual transport and optical properties suggesting that they may have a 
gap at the Fermi leveLtOil^ This observation establishes a possible new class of semiconduct- 
ing intermetallics, which we have identified as good subjects for the theoretical study of band 
gap formation in intermetallics and its relation to structural stability.tJ These compounds 
crystallize in the MgAgAs structure and are related to the metallic Heusler compounds 
Ni2SnM by the removal of an fee Ni sublattice.ll3 This simple crystal structure makes them 
quite suitable for a first principles approachta to the microscopic understanding of their un- 
usual optical and electronic properties. While a complete understanding of these properties 
would require incorporating the effects of a high concentration of defects, calculations for 
the ideal structures should provide a useful starting point for further analysis. To under- 
stand the origin of the semiconducting gap in this system, one needs more than "simple" 
arguments. From a tight-binding analysis of the first principles bands and charge densities, 
both the pd hybridization on the SnM rocksalt substructure and the symmetry breaking 
due to the Ni sublattice are found to be crucial factors in gapping the electronic spectrum. 
Examination of total energies of ternary compounds and hypothetical binary substructures, 
calculated from first principles, allows us to conclude that the formation of the gap is not 
directly associated with structural stability in this system. Finally, the use of first princi- 
ples information also assists us in identifying chemical trends for the prediction of related 
semiconducting intermetallics. 

The rest of the paper is organized as follows. In Section II, we briefly discuss the technical 
details of the ab initio pseudopotential calculations for the NiSnM and Ni2SnM compounds. 



In Section III, we first examine their ground state structural properties and identify the 
correct atomic arrangement for the NiSnM compounds. This is followed by the results of 
band-structure and DOS calculations which characterize, within the local density approx- 
imation (LDA), the Ni2SnM compounds as normal metals, and the NiSnM compounds as 
narrow gap semiconductors with indirect gaps near 0.5 eV. In Section IV, we discuss these 
results in a tight-binding framework and by considering decompositions of the ternary struc- 
tures into binary substructures plus a third element sublattice. Particular attention is given 
to the origin of the gap and structural stability of NiSnZr. We compare our results with 
available experimental measurements and, using the lessons from our analysis, we predict a 
new semiconducting intermetallic compound. Finally, our results are summarized in Section 
V. 

II. CALCULATIONAL PROCEDURE 

For the calculation of band structures and total energies, we used the ab initio pseu- 
dopotential method with a plane wave basis set and the conjugate gradients algorithm.li3 
The general technical issues that arise in the application of this method to intermetallic 
systems containing transition metals have been discussed in detail in Ref. [l^. In this study, 
for Sn, Zr and Hf, we used the scalar- relativistic pseudopotentials of Bachelet, Hamann, 
and Schliiterll2l (BHS). For Ti and Ni, we constructed Hamann, Schliiter, and Chiang (HSC) 
pseudopotentialsEy for the s and p orbitals, and optimized pseudopotentialsB (OPT) for the 
d orbitals. The exchange-correlation potential used in these calculations is the functional of 
Ceperley and Alder as parametrized by Perdew and Zunger.E3 Complete information con- 
cerning pseudopotential constuction is given in Table I. These pseudopotentials were put 
into separable form using single or double projectors for each angular momentum (Table I) 
with the / = component as the local potential.c3 In all cases, the absence of ghost states 
and the transferability of pseudopotentials were checked by applying the ghost theorem of 
Gonze et alB-i and comparing the logarithmic derivatives for these nonlocal pseudopotentials 



with their all-electron values over a wide energy range. 

All total energy and band-structure calculations were done on IBM RS/6000 workstations 
using the conjugate gradients program castep 2.1,1^ slightly modified for metallic systems. 
The plane wave cutoff used in the calculations was determined by requiring the kinetic 
energy in the Fourier components of all atomic pseudo wave functions to be less than 1 
mRy. This cutoff was set by the sharply peaked 3d valence states of Ni to be 52 Ry, allowing 
us to perform accurate calculations with 2200-2800 plane waves depending on the size of the 
unit cell. Due to the metallic nature of some of the compounds investigated in this work, 
we used a Fermi-Dirac broadeningc3 of 50 meV for k-point sampling with Monkhorst-PackE3 
grids of g = 8 and q = 10, resulting in 10 and 19 k points in the irreducible Brillouin zone 
(BZ), respectively. Lattice constants and bulk moduli were obtained by fitting to the Birch 
form, as described in Ref. p!8| . 

The band-structure calculations for the NiSnM and Ni2SnM compounds were performed 
with self-consistent charge densities computed with a Monkhorst-Pack grid oi q = 10. For 
DOS calculations, we first used this charge density to evaluate the band structure on a 
Monkhorst-Pack grid oi q = 20 (110 k points in the irreducible BZ), and interpolated the 
eigenvalues to a grid of ~ 225 000 k points in the irreducible fee BZ using the interpolation 
scheme of Monkhorst and PackBZl based on a global BZ fit. We then applied the Gilat- 
Raubenheimer methodc3 to this fine mesh in the irreducible BZ. To check the accuracy of 
this interpolation, the calculated band structures along certain high symmetry lines were 
compared with the interpolated values, and the agreement was found to be very good within 
a few percent. Finally, for the band-structure plots and for use in the tight-binding analysis, 
symmetry labels along high symmetry lines were assigned using projection operators for the 
corresponding irreducible representations. 



III. RESULTS 

The ternary crystal structures included in this study are the BiFs structure for the 
Ni2SnM compounds and the MgAgAs structure for the NiSnM conipounds.E3 The crystal 
structure of Ni2SnM (space group FmSm), shown in Fig. 1(a), can be viewed as a rocksalt 
arrangement of Sn and M atoms with Ni atoms at the center of each cube of Sn-M atoms. 
The noncentrosymmetric crystal structure of NiSnM (space group FASm), shown in Fig. 
1(b), can be obtained from that of Ni2SnM by eliminating one of the fee Ni sublattices, 
resulting in tetrahedral coordination of the Sn and M atoms by Ni. We also considered two 
additional distinct MgAgAs type crystal structures obtained by interchanging the Ni, Sn and 
M sublattices. These three different crystal structures can be distinguished by identifying 
the doubly tetrahedrally coordinated element; the relevant coordination environments are 
shown in Fig. 2(a). We refer to these three structures as the a, (3, and 7 phases depending 
on whether Ni, Sn, or M is doubly tetrahedrally coordinated, respectively.E3 

First, we performed total energy calculations at several lattice constants for all six com- 
pounds Ni2SnM and NiSnM (in the a phase) to determine their ground state structural 
properties. The results are summarized in Table II. As seen from this table, the calculated 
equilibrium lattice constants are in very good agreement with experiment, with slight un- 
derestimates typical of the LDA. To our knowledge, no measurements of the bulk moduli for 
these compounds have been reported in the literature; the values in the table are therefore 
predictions for the bulk moduli of the Ni2SnM and NiSnM compounds. We also performed 
total energy calculations for the a, /3, and 7 phases of NiSnM. Figure 2(b) shows the results 
of these calculations for M = Zr, and the a phase is seen to be the lowest energy struc- 
ture (2.28 and 2.85 eV/unit cell below the (3 and 7 phases, respectively), consistent with 
reported structural determinations. The experimental lattice constant for the NiSnZr com- 
pound (6.11 A)Ej is also closest to the LDA lattice constant of the a phase. This calculation 
therefore unambiguously determines the crystal structure of the NiSnM compounds as the a 
phase. The reasons behind the relative stability of the ct, /3, and 7 phases will be discussed 



in Section IV. 

Next, we performed band-structure and DOS calculations for both Ni2SnM and NiSnM 
compounds as a first step towards a microscopic understanding of their electronic structures. 
Figure 3 shows the calculated band structure along high symmetry lines of the fee BZ for 
Ni2SnZr. This exhibits normal metallic behavior with bands crossing the Fermi level along 
various directions, which results in a finite DOS at the Fermi level, as shown in Fig. 4. Very 
similar band structures and DOS were also obtained for Ni2SnTi and Ni2SnHf. On the other 
hand, band-structure calculations for the lowest energy structure of NiSnM (the a phase) 
show that the Fermi level does not intersect any bands along the high symmetry directions, 
as can be seen in Fig. 5 for the case of NiSnZr. In fact, DOS calculations show that for these 
compounds there is a gap in the electronic DOS at the Fermi level (Fig. 6). For all three 
compounds, this is an indirect Y -^ X gap near 0.5 eV, while the direct gap, close to 1.1 
eV, occurs at X (Table III). In contrast, band-structure calculations for the higher energy /? 
and 7 phases of NiSnM show metallic behavior (Fig. 7). Therefore, these LDA calculations 
characterize the NiSnM compounds in the a phase as narrow gap semiconductors, while 
Ni2SnM compounds and the /? and 7 phases of NiSnM are found to be normal metals. 

IV. DISCUSSION 

Using the results of first principles total energy and band-structure calculations, we can 
develop a physical understanding of the microscopic mechanisms in the formation of the 
band gap in the a phase of NiSnM, and of the differences between these compounds and 
related systems, specifically the /3 and 7 phases and Ni2SnM. In Section IV A, we investigate 
the origin of the band gap in a tight-binding framework. In Section IV B, we consider 
the three phases of the ternary compounds NiSnM using the decompositions into a binary 
substructure plus a third element sublattice, and show that these decompositions can be 
quite useful in discussing the stability and electronic structures of ternary compounds. After 
a brief discussion of the differences between the metallic Heusler compounds Ni2SnM and the 
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semiconducting NiSnM compounds in Section IV C, we examine the available experimental 
data on the NiSnM compounds from the point of view of our theoretical calculations in 
Section IV D. The relationships between band gap formation and structural stability in 
intermetallic compounds are considered in Section IV E, and the applicability of these ideas 
to the NiSnM compounds is examined. Finally, in Section IV F, we use the lessons gained 
from this analysis to predict a new semiconducting intermetallic compound. We note that 
for definiteness, the discussion in these sections is for M = Zr, though the overall features 
are the same for M = Ti and Hf. 

A. Origin of the band gap 

Our tight-binding analysis begins with identification of the orbital characters of the va- 
lence and conduction bands. With a plane- wave basis set, this is accomplished by examining 
contour plots of the charge densities of individual states at various k points. We first consider 
the bands well below the Fermi level. The lowest-lying band in Fig. 5 has predominantly 
Sn s orbital character, while the narrow bands above it have predominantly Ni d character. 
This seems to suggest a picture in which bands arising from the Sn s, Ni d, and either Sn 
p or t2g type Zr d orbitals are fully occupied, forming a closed shell system with a gap to 
the next excited state. We therefore examined the triply degenerate Fis and the doubly 
degenerate X5 states just below the Fermi level, shown in Fig. 8(a). The charge density 
plots in the (110) plane show that the Fis states have mainly Zr d character [Fig. 8(c)], 
while the X^ states have mainly Sn p character [Fig. 8(e)]. This change in orbital character 
shows that the gap in NiSnM compounds cannot be due to a simple closed shell picture, 
since in that case we would expect the highest occupied band to be derived from either Sn 
p or Zr d orbitals. 

A detailed analysis of band behavior is therefore needed to explain the formation of 
the gap. We focus on the bands F — A — A as capturing the essential features of the 
gap mechanisms. Symmetry labels and orbital characters allow the separation of bands 



into different subspaces which can be considered one at a time. We start with the highest 
occupied valence band A3 4. The doubly degenerate A3 4 conduction band from Fis to X^ 
[Fig. 8(a)] has the same symmetry label, but complementary orbital character, changing 
from mainly Sn jo at F [Fig. 8(b)] to mainly 7^i d aX X [Fig. 8(d)]. To describe this pair of 
bands, we use a tight-binding basis including the Sn p^, py, and Zr d^z, dy^ orbitals.E3 The 
behavior of these bands, shown with thick lines in Fig. 8(a), can be explained as arising from 
dispersive Zr d and Sn p type bands which hybridize along A and anticross. The reason that 
this does not look like a familiar anticrossing of two bands (where the bands approach and 
then slightly repel each other) is that this pd hybridization is large and strongly k-dependent, 
with the maximum value halfway along A. Specifically, starting with nonzero interactions 
only for Sn-Sn and Zr-Zr nearest neighbors in the fee lattice, the dispersion Ex^y{a) for the 
Sn Pj. and Py orbitals is 

Ex,y{ct) = €p + 2{ppa + pp-Ti) + {2ppa + Qppn) cos 7ra (4.1) 

while the dispersion Exz,yz{o!) for the Zr d^z and dyz orbitals is 

Exz,yz{.Oi) = td + 2{dd'K + dd6) + {3dda + 2ddn + 3dd6) costto (4.2) 

where a runs from (at F) to 1 (at X). In these equations, e^ and e^ denote the onsite 
orbital energies of the p and d orbitals, respectively, and the definitions of the tight-binding 
parameters are as in Ref. ^. If the parameters are taken to have the usual approximate 
relative magnitudes and signs {ddrc ^ —^dda, ddS ~ 0, ddcr < 0, pp-n ^ —\ppo^, and 
ppa > 0)f^ the p type bands disperse downward with a maximum at F, while the d type 
bands disperse upward with a minimum at F. With a suitable choice of the parameters, the 
two bands will cross. However, including a nearest neighbor Sn-Zr pd interaction will result 
in an anticrossing of these two bands, with a matrix element 2pdn sin ira for this geometry, 
which is largest halfway along A. This tight-binding parametrization reproduces the first 
principles dispersion and is consistent with the orbital character of the states obtained from 
charge density plots. This shows that the strong pd hybridization in the SnZr rocksalt 
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substructure is a necessary ingredient in the formation of the gap, for if it were zero, the ri5 
states below (above) the Fermi level would join with the X5 states above (below) the Fermi 
level resulting in a finite DOS at the Fermi level. However, this hyridization is not in itself 
sufficient for the gap formation in the ternary NiSnZr compound. 

To complete the analysis of gap formation in the a phase of NiSnZr, we must consider 
the conduction band minimum X3 [Fig. 8(a)], which has mainly Zr d^y character. Since the 
valence band maximum Fis also has Zr d^y character, an unhybridized Zr dxy band, with Ai 
symmetry, must anticross with Ai bands derived from other orbitals for the gap to form. 
However, this anticrossing cannot be explained by considering only the SnZr substructure 
as in the case of A3 4 bands. Due to the Oh symmetry of the rocksalt substructure, the 
matrix elements between the Zr d^y band along A and the bands derived from all other Zr d 
and Sn p orbitals vanish, resulting in zero anticrossing. For the gap to form, it is necessary 
to include the Ni d orbitals, lowering the symmetry of the system to T^. Although the Ni 
d states lie well below the Fermi level, the Ni d^y and d^-^2_^2 orbitals mediate an indirect 
interaction between the Zr d^y band and the bands derived from the Zr d^z2_j.2 and Sn p^ 
orbitals. These interactions lead to a nonzero anticrossing of these bands, opening the gap 
at the Fermi level. Therefore, the symmetry breaking produced by the third element (Ni) 
sublattice is the essential additional ingredient in the formation of the gap. 

B. Ternary compounds from binary substructures 

The analysis of the previous section showed that the states near the Fermi level in the 
a phase of NiSnZr are dominated by the atomic orbitals of the SnZr rocksalt substructure. 
This suggests an approach to the analysis of the electronic structure and stability of the 
ternary compound through the identification of a key binary substructure. In this approach, 
the essential features of electronic structure and stability should already be evident in the 
binary substructure, while addition of the third element sublattice leads to recognizable 
modifications such as (i) introduction of extra orbitals hybridizing with the orbitals of the 
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binary substructure, (ii) change in the symmetry of the crystal structure, and (iii) charge 
transfer between the binary substructure and the third element sublattice, and other changes 
in the self-consistent charge density. 

We apply this approach to NiSnZr by computing the total energies and band structures 
for the binary substructures appearing in the decomposition of the a, /3 and 7 phases. There 
are six distinct substructures: SnZr, NiZr and NiSn in the rocksalt (-B1) and zincblende (-B3) 
structures. The total energies of these six substructures at the lattice constant of the a phase, 
with lists of the corresponding NiSnZr phases, are given in Table IV. For each of the three 
binary compositions, the lower energy structure (between Bl and B3) is the one contained 
in the decomposition of the a phase. However, consideration of the magnitudes of the BI- 
BS energy differences singles out the decomposition containing the SnZr substructure. The 
Bl structure of SnZr is 2.03 eV lower in energy than the B3 structure. This value is quite 
close to the energy difference between the a phase, containing the Bl SnZr substructure, 
and the /3 and 7 phases, which contain the B3 SnZr substructure (2.28 eV and 2.85 eV, 
respectively). In contrast, the B1-B3 energy differences for NiZr and NiSn are too small 
(0.26 eV and 0.23 eV per unit cell, respectively) to be relevant to the stabilization of the a 
phase relative to the (3 and 7 phases. Therefore, the low energy of the NiSnZr compound in 
the a phase can be mainly attributed to the stability of the SnZr Bl substructure. 

There also appear to be useful relationships between the band structures of the ternary 
compounds and those of binary substructures. For example, the a phase of NiSnZr, where 
the SnZr substructure is rocksalt, has a gap while the (3 and 7 phases, with SnZr in a 
zincblende arrangement, are both metallic. We investigate these relationships in more detail 
by examining the band structures of the six substructures (Fig. 9) and comparing them with 
the band structures of the ternary a [Fig. 8(a)], (3 [Fig. 7(a)] and 7 [Fig. 7(b)] phases. In 
the B3 phases of NiSn and NiZr, the crystal field splittings of the Ni and Zr d states at F are 
noticeably larger than in the Bl phases. The reason is that the point symmetry is lowered 
from Oh in the Bl structure to T^ in the B3 structure, resulting in a mutual repulsion of t2g 
type d and p orbitals in the latter. Comparison of the Bl with the B3 band structures for 
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the three compositions shows a striking difference between NiZr and NiSn, on the one hand, 
and SnZr, on the other hand, in the position of the Fermi level and dispersion of the nearby 
bands. The very small change in the position of the Fermi levels for the two structures 
of NiZr and NiSn is not surprising, because the Ni d orbitals are rather low-lying, and do 
not substantially affect states around the Fermi level in these compounds. In contrast, the 
arrangement of the bands and position of the Fermi level changes noticeably from SnZr Bl 
to SnZr B3. This occurs because Zr d and Sn p orbitals are partially occupied orbitals 
involved in bonding, making the band dispersions and the positions of the Fermi level quite 
sensitive to the atomic arrangement. 

We can also compare the band structures of the three phases of the ternary compound 
with the band structures of the Bl and B3 structures of SnZr that appear in their decom- 
positions. For example, comparison of Fig. 8(a) and the band structure of SnZr in the Bl 
structure (Fig. 9) shows that the low-lying Sn s band, the conduction bands well above 
the Fermi level, and the bands arising from the Sn p^, py and Zr d^z, dyz hybridizations 
(labelled by A5 in Fig. 9) look quite similar to those in Fig. 8(a). On the other hand, the 
Ai band in Fig. 9, arising from the doubly degenerate ri2 state and ending up as the X2 
state at X is substantially modified by the introduction of the Ni sublattice. In the ternary 
compound, this band anticrosses due to the presence of Ni d orbitals and reduced point 
group symmetry. Closer examination shows that the introduction of the Ni d orbital greatly 
increases the crystal field splitting of the Zr d states at F, with the Fermi level in the a phase 
of NiSnZr placed just inside this splitting. Thus, within the substructure decomposition ap- 
proach, we can account for the formation of the gap. The band structure of NiSnZr in the 
7 phase [Fig. 7(b)] also has similarities with the band structure of SnZr in the B3 structure 
(Fig. 9). Since the Sn p— Zr d interaction pushes the t2g type Zr d states down in energy 
in the B3 structure, the introduction of Ni d orbitals results in the formation of a d-band 
complex below the Fermi level, while the downward dispersing bands well above the Fermi 
level are not substantially affected by this. On the other hand, the relationship between the 
band structure of the (3 phase of NiSnZr [Fig. 7(a)] with the band structure of its B3 SnZr 
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substructure is much more complicated, because introducing the Ni d orbitals changes the 
mixing of the p and d orbitals to a greater extent due to double tetrahedral coordination of 
Sn. In this case, it is not possible to treat the effects of the Ni d orbitals perturbatively. 

C. Comparison of Ni2SnM with NiSnM 

The Ni2SnM compounds can be produced from the a phase of NiSnM by adding a Ni 
sublattice, which changes the system from a narrow gap semiconductor to a metal. The effect 
of the additional sublattice responsible for this change in electronic structure is the raising 
of the noncentrosymmetric Td symmetry to the centrosymmetric Oh point group symmetry, 
with consequent changes in the symmetry properties of atomic-orbital Bloch functions. In 
particular, if (p denotes any of the five d or s orbitals of Ni, the functions that transform 
according to the irreducible representations of the point group of any high symmetry point 
in the BZ have the form -^((p^^^ =p <P^'^^), where 0*^*-* (z = 1,2) denotes the (p orbital of the 
i^^ Ni atom in the unit cell. Also, any orbital that transforms according to Ai (X3) in the 
noncentrosymmetric compound transforms in Ni2SnM according to either Ai or Ag (X4 or 
X2). Therefore, the Ai and A2 bands near the Fermi level (Fig. 3) do not anticross as 
they do in NiSnM, where they both carry the Ai label. As a result, the Ai band cuts the 
Fermi level, producing the metallic behavior. Finally, in our calculations we do not see any 
evidence of a vacancy banda proposed to explain the experimental data for NiSnM. This 
is not too surprising for the perfectly clean samples on which our calculations were based, 
because unlike the case of heavy doped semiconductors, the vacancies in this system are 
very well ordered (in an fee lattice) and have a high concentration. 

D. Comparison with experiment 

What light do these calculations shed on the experimentally observed electronic and 
optical properties of the NiSnM compounds? First of all, the very high low-temperature 
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resistivities with a simple exponential temperature dependence observed in the NiSnM com- 
pounds are consistent with the the appearance of a semiconducting gap in the calculated 
band structures. However, specific heat measurements show a small but finite DOS at the 
Fermi level.li3 Furthermore, the calculated gaps are all near 0.5 eV. These are more than 
twice as large as the gap parameters 0.12, 0.19, and 0.22 eV for Ti-, Zr-, and Hf-based 
compounds, respectively, extracted from resistivity measurementsEj and features in the op- 
tical refiectivity R{io)f-^ While LDA errors (either underestimates or overestimates) in the 
values of the band gaps are expected, a more serious problem which could give rise to these 
discrepancies is poor sample quality, particularly substitutional disorder. In NiSnZr, there is 
direct evidence from x-ray diffraction for Sn-Zr substitutional disorder at the 10-30% level.ll3 
The degree of substitutional disorder seems to be correlated with changes in the electronic 
properties. One example of this is that the DOS at the Fermi level appears to decrease as the 
degree of disorder decreases. In particular, the zero temperature resistivity increases dra- 
matically upon annealing the samples.t^ Also, the DOS at the Fermi level from specific heat 
measurement st3 is lowest for the Ti-based compound, which has the smallest disorder in the 
position of the Ti and Sn atoms, indicating that a complete gap may be expected to open in 
a very clean sample. To improve our understanding of the effect of the Sn-Zr substitutional 
disorder on the size of the gap, we performed total energy and band-structure calculations 
for Ni(Sni_^Zr^.)(Sna.Zri_^.) for a few values of x between and 0.5 in the virtual crystal 
approximation. At around x = 0.15, the gap closes to form a semimetal, with the valence 
band maximum (conduction band minimum) moving up (down) in energy and touching the 
Fermi level. Therefore, to achieve experimental gap values comparable to those calculated 
for the ideal crystals, special attention to sample preparation is needed. In such well-ordered 
samples, reflectivity measurements would have to be extended to higher frequencies than 
those used in previous studies (in Ref. ITl| cj < 5000 cm'^ ) to see the onset of absorption. 
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E. Band gap and structural stability 

Direct relations between the formation of a (pseudo)gap and the stabihty of the structure 
have been proposed in several contexts, including the Hume Rothery rule,0'Q Jahn- Teller 
effects in molecules,Eil and the "coloring" patterns of atoms over the sites of a given latticeS 
If the formation of a gap in the DOS is an important factor in the structural stability of 
intermetallic compounds, the expectation is that semiconducting intermetallics would be 
fairly common, while, as we have already discussed, they are quite rare. Therefore, it is 
useful to examine this issue in the present case. 

For NiSnM compounds in their a, (3, and 7 phases, it is true that the structure with 
the gap has the lowest energy, and there is a correlation between the DOS at the Fermi 
level and the structural energies of the three phases. However, the present calculations 
suggest that the formation of the gap is not directly responsible for the structural stability. 
Instead, the stability can be understood from the substructure decomposition approach. 
As discussed in Sec. IV B, the energy of the SnZr Bl structure is 2.03 eV per unit cell 
lower than that of the B3 structure. From Fig. 9, it is clear that neither the Bl nor the 
B3 SnZr structures have a gap. The ternary phases are obtained by incorporating the Ni 
sublattice into the binary SnZr substructures. This incorporation produces a gap in the a 
phase, as a result of symmetry breaking, but leads to only a small change in the relative 
energy (from 2.03 eV per unit cell to energy differences of 2.28 eV and 2.85 eV per unit 
cell between the a and the (3 and 7 phases, respectively). Therefore, the stabilization of 
the lowest energy phase is primarily due to the relative stability of the binary SnZr Bl 
substructure, while the contribution to stabilization associated with the incorporation of 
the third element (Ni) sublattice, which results in the formation of the gap, is relatively 
minor. This conclusion is not inconsistent with the relations between gap formation and 
stability observed elsewhere. In the Hume-Rothery rule, a gap in the DOS establishes 
the limit of stability of a particular alloy structure as a function of electron count, rather 
than an enhanced stability specifically of the semiconducting phase. In the Jahn- Teller and 
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coloring analyses, the correlation between structural stability and gap formation is identified 
in classes of structures which are otherwise extremely similar (related by small structural 
relaxations in the Jahn- Teller case and possessing identical nearest neighbor environments 
for the cations in the coloring problem). However, in the NiSnM compounds, the a, /? 
and 7 phases are related by sublattice interchanges, not small displacements, resulting in 
different nearest neighbor environments in each phase. It is therefore reasonable that the 
substructure decomposition approach, with the identification of the low-energy rocksalt SnZr 
substructure, should be a more useful indicator of structural stability than gap formation 
for these compounds. 

F. Considerations for ne^v semiconducting intermetallics 

From these results, the simple chemical and structural trends in gap formation can be 
identified and used to assist in the prediction of new semiconducting intermetallics. We have 
already seen that a closed shell rule with 18 electrons is too simple. The ternary nature of 
the compounds, particularly the electronegativity difference between the constituents of the 
rocksalt substructure, and the structural arrangement, particularly the double tetrahedral 
coordination of the late transition metal, seem to be important factors, since the hypothetical 
compounds NiSn2 and NiZr2, and the /3 and 7 phases of NiSnM compounds are metallicE^ 
On the other hand, the empirical rule that the total number of d electrons from the two 
transition metals in their ground states should be 10 works for the NiSnM compounds as 
well as for the compounds formed by the column substitutions of Ni with Pd and Pt.ll3 It is 
also necessary that the transition metal that is doubly tetrahedrally coordinated should have 
a significantly larger number of d electrons in the atomic ground state, so that the bands 
originating from these d orbitals will lie deep below the Fermi level, and not destroy the 
pd hybridizations by mixing significantly with the d orbitals of the other transition metal. 
CoVSn in the MgAgAs structure with Co doubly tetrahedrally coordinated by Sn and V 
satisfies these rules. First-principles band structure calculations confirm that CoVSn has an 
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L ^ X gap of 0.8 eV at the Fermi level (Fig. 10). This compound is not reported to exist in 
the intermetallic compound database. However, there is a good chance that this might be a 
stable compound based on the analogy with the a phase of NiSnM compounds. Therefore, 
we predict that if this compound can successfully be prepared in this structure, it will be a 
semiconducting intermetallic. 

V. CONCLUSIONS 

In this paper, we used the pseudopotential total energy method to investigate the stability 
and electronic structure of the ternary intermetallic compounds NiSnM (M = Ti, Zr, Hf). 
Our calculations revealed an indirect semiconducting gap near 0.5 eV for all three compounds 
in the lowest energy a phase, while the inversion symmetric Ni2SnM compounds and the 
higher energy (3 and 7 phases of NiSnM were found to be metallic. From a tight-binding 
analysis of the first-principles bands and charge densities, both the pd hybridization on 
the SnZr rocksalt substructure and the symmetry breaking due to the Ni sublattice were 
found to be crucial factors in gapping the electronic spectrum of a-NiSnZr. By considering 
all binary substructures which appear in the decomposition of the three phases of NiSnZr, 
we concluded that the stabilization of the lowest energy a phase is primarily due to the 
relative stability of the binary SnZr Bl substructure, while the contribution to stabilization 
associated with the incorporation of the third element (Ni) sublattice, which results in the 
formation of the gap, is relatively minor. Analysis of the band structures for the narrow 
gap semiconductor NiSnZr and the metallic Ni2SnZr showed that the raising of the point 
group symmetry from the noncentrosymmetric T^ to the centrosymmetric Oh group by the 
addition of the Ni sublattice is responsible for this difference in the electronic properties. 
Our theoretical findings are consistent with experimental observations, if poor sample quality 
and substitutional disorder in the samples are taken into account. Finally, we predicted that 
the ternary compound CoVSn in the MgAgAs structure with an SnV rocksalt substructure 
will form a semiconducting intermetallic, if it can be prepared successfully in this structure. 
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FIGURES 

FIG. 1. The crystal structures of (a) Ni2SnM; (b) NiSnM compounds. 

FIG. 2. (a) Three possible coordination environments corresponding to the a (bottom), (3 
(middle), and 7 (top) phases of NiSnM; (b) Total energy versus the lattice constant for the a, [3, 
and 7 phases of NiSnZr. The energies shown in this graph are with respect to the minimum energy 
corresponding to the LDA lattice constant of the a phase. 

FIG. 3. The energy bands of Ni2SnZr along the high symmetry lines of the fee Brillouin zone. 
The calculations were performed at the experimental lattice constant of 6.27 A. The dotted line is 
the Fermi level. The symmetry labels are shown only for states that are discussed in Section IV C. 

FIG. 4. Electronic density of states for Ni2SnZr. The calculations were performed at the 
experimental lattice constant of 6.27 A. The dotted line is the the Fermi level. 

FIG. 5. The energy bands of NiSnZr along the high symmetry lines of the fee Brillouin zone. 
The calculations were performed at the experimental lattice constant of 6.11 A. The dotted line is 
the Fermi level. 

FIG. 6. Electronic density of states for NiSnZr. The calculations were performed at the exper- 
imental lattice constant of 6.11 A. The dotted line is the the Fermi level. 

FIG. 7. The energy bands of NiSnZr in the (a) /3; (b) 7 phase along A. The calculations were 
performed at a lattice constant of 6.11 A. The dotted lines represent the position of the Fermi 
level. 

FIG. 8. (a) The energy bands of NiSnZr along A. The thick curves are the doubly degenerate 
(due to time reversal symmetry) A3^4 bands mentioned in the text; (b)-(c) Charge density contour 
plots for the triply degenerate Fis states; (d)-(e) Charge density contour plots for the doubly 
degenerate X^ states. The circles filled with x marks and thick dots represent the positions of the 
Sn and Zr atoms, respectively. 
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FIG. 9. The energy bands for the SnZr, NiZr, and NiSn binary compounds in the Bl and B3 
structures along A. All calculations were performed at a fixed lattice constant of 6.11 A. The 
dotted lines represent the position of the Fermi level. 

FIG. 10. Band structure of CoVSn in the MgAgAs structure along L — A — T — A — X. The 
dotted line is the Fermi level. The calculations were performed at a lattice constant of 5.9 A. 
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TABLES 
TABLE L Parameters and methods used in the construction of the pseudopotentials. For OPT, 

qc and N in the last two columns refer to the cutoff above which the kinetic energy is minimized 

and the number of Bessel functions used to expand the pseudo wave functions, respectively.^^ 

Element Configuration Type Projector^ r^ (a.u) q^, (Ry) N 

Ni 4s°-^5 4p0-25 Sd^ s, p HSC p, single 1.2, 1.2 

d, OPT d, single 2.0 7.15 10 

Ti 4s0-^5 4p0-25 3d2 s, p HSC p, single 1.3, 1.4 

d, OPT d, single 1.9 6.80 4 

Sn 5si 5p0-5 5d0-5 BHS p, double 

d, double 

Zr SsO'^^ 5p°-25 Ad? BHS p, double 

d, double 

Hf esO'^^ 6p°-25 5d2 BHS p, single 

d,double 

'^The I = component was chosen as the local pseudopotential. 
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TABLE II. The experimental and calculated structural properties of the Ni2SnM and NiSnM 
compounds. 

Compound &exp. (A) ^ a-ca/c. (^) B (Mbar)*^ 

NisSnTi 6.09 6.00 1.64 (2.14) 

Ni2SnZr 6.27 6.24 1.49 (1.59) 

NisSnHf 6.24 6.13 1.36 (1.72) 

NiSnTi 5.92 5.814 1.11 (1.43) 

NiSnZr 6.11 6.07 1.19 (1.29) 

NiSnHf 6.066 5.95 1.10 (1.41) 
^Prom Ref. |l| and Fig. 5 in Ref. |l|. 
''The values in parentheses are the bulk moduli calculated at the LDA lattice constant. 

TABLE III. LDA band gaps for NiSnM compounds. The indirect gap occurs along T to X, 
while the direct gap occurs at X. 

Compound Indirect gap (eV) Direct gap (eV) 

NiSnTi 0.51 1.20 

NiSnZr 0.51 1.05 

NiSnHf 0.48 1.16 
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TABLE IV. Total energies for six different binary compounds which appear in the decompo- 
sition of the a, /3, and 7 phases of NiSnZr into Bl and B3 structures. The Greek letters in the 
parentheses under the second column give the phase of NiSnZr which will be obtained by the ad- 
dition of the third type of element into either the tetrahedral or the octahedral hole of the fee unit 
cell. All calculations were performed at a fixed lattice constant of 6.11 A. 

Compound Structure Total Energy (eV) 

SnZr Bl (a) -180.200 

SnZr B3 (/3,7) -178.172 

NiZr Bl (/?) -997.115 

NiZr B3 (0,7) -997.375 

NiSn Bl (7) -1011.445 

NiSn B3 (a,d) -1011.671 
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